Pediatric Chordoma: A Tale of Two Genomes.

Little is known regarding the genomic alterations in chordoma, with the exception of loss of SMARCB1, a core member of the SWI/SNF complex, in poorly differentiated chordomas. A TBXT duplication and rs2305089 polymorphism, located at 6q27, are known genetic susceptibility loci. A comprehensive genomic analysis of the nuclear and mitochondrial genomes in pediatric chordoma has not yet been reported. In this study, we performed whole exome and mitochondrial DNA (mtDNA) genome sequencing on 29 chordomas from 23 pediatric patients. Findings were compared with that from whole genome sequencing datasets of 80 adult skull base chordoma patients. In the pediatric chordoma cohort, 81% percent of the somatic mtDNA mutations were observed in NADH complex genes, which is significantly enriched compared to the rest of the mtDNA genes (p=0.001). In adult chordomas, mtDNA mutations were also enriched in the NADH complex genes (p<0.0001). Furthermore, a progressive increase in heteroplasmy of non-synonymous mtDNA mutations was noted in patients with multiple tumors (p=0.0007). In the nuclear genome, rare likely germline in-frame indels in ARID1B, a member of the SWI/SNF complex located at 6q25.3, were observed in five pediatric patients (22%) and four patients in the adult cohort (5%). The frequency of rare ARID1B indels in the pediatric cohort is significantly higher than that of the adult cohort (p=0.0236, Fisher's exact test), but they were both significantly higher than that in the ethnicity-matched populations (p<5.9e-07 and p<0.0001174, respectively). Implications: germline ARID1B indels and mtDNA aberrations appear important for chordoma genesis, especially in pediatric chordoma.

An Exome Capture-Based RNA-Sequencing Assay for Genome-Wide Identification and Prioritization of Clinically Important Fusions in Pediatric Tumors.

This study reports the development of an exome capture-based RNA-sequencing assay to detect recurring and novel fusions in hematologic, solid, and central nervous system tumors. The assay used Twist Comprehensive Exome capture with either fresh or formalin-fixed samples and a bioinformatic platform that provides fusion detection, prioritization, and downstream curation. A minimum of 50 million uniquely mapped reads, a consensus read alignment/fusion calling approach using four callers (Arriba, FusionCatcher, STAR-Fusion, and Dragen), and custom software were used to integrate, annotate, and rank the candidate fusion calls. In an evaluation of 50 samples, the number of calls varied substantially by caller, from a mean of 24.8 with STAR-Fusion to 259.6 with FusionCatcher; only 1.1% of calls were made by all four callers. Therefore a filtering and ranking algorithm was developed based on multiple criteria, including number of supporting reads, calling consensus, genes involved, and cross-reference against databases of known cancer-associated or likely false-positive fusions. This approach was highly effective in pinpointing known clinically relevant fusions, ranking them first in 47 of 50 samples (94%). Detection of pathogenic gene fusions in three diagnostically challenging cases highlights the importance of a genome-wide and nontargeted method for fusion detection in pediatric cancer.

Low-pass whole-genome and targeted sequencing of cell-free DNA from cerebrospinal fluid in pediatric patients with central nervous system tumors.

Background: Central nervous system tumors are the most common pediatric solid tumors and the most frequent cause of cancer-related morbidity in childhood. Significant advances in understanding the molecular features of these tumors have facilitated the development of liquid biopsy assays that may aid in diagnosis and monitoring response to therapy. In this report, we describe our comprehensive liquid biopsy platform for detection of genome-wide copy number aberrations, sequence variants, and gene fusions using cerebrospinal fluid (CSF) from pediatric patients with brain, spinal cord, and peripheral nervous system tumors. Methods: Cell-free DNA was isolated from the CSF from 55 patients, including 47 patients with tumors and 8 controls. Results: Abnormalities in cell-free DNA were detected in 24 (51%) patients including 11 with copy number alterations, 9 with sequence variants, and 7 with KIAA1549::BRAF fusions. Positive findings were obtained in patients spanning histologic subtypes, tumor grades, and anatomic locations. Conclusions: This study demonstrates the feasibility of employing this platform in routine clinical care in upfront diagnostic and monitoring settings. Future studies are required to determine the utility of this approach for assessing response to therapy and long-term surveillance.

Resilience to autosomal dominant Alzheimer's disease in a Reelin-COLBOS heterozygous man.

We characterized the world's second case with ascertained extreme resilience to autosomal dominant Alzheimer's disease (ADAD). Side-by-side comparisons of this male case and the previously reported female case with ADAD homozygote for the APOE3 Christchurch (APOECh) variant allowed us to discern common features. The male remained cognitively intact until 67 years of age despite carrying a PSEN1-E280A mutation. Like the APOECh carrier, he had extremely elevated amyloid plaque burden and limited entorhinal Tau tangle burden. He did not carry the APOECh variant but was heterozygous for a rare variant in RELN (H3447R, termed COLBOS after the Colombia-Boston biomarker research study), a ligand that like apolipoprotein E binds to the VLDLr and APOEr2 receptors. RELN-COLBOS is a gain-of-function variant showing stronger ability to activate its canonical protein target Dab1 and reduce human Tau phosphorylation in a knockin mouse. A genetic variant in a case protected from ADAD suggests a role for RELN signaling in resilience to dementia.

Combined low-pass whole genome and targeted sequencing in liquid biopsies for pediatric solid tumors.

We designed a liquid biopsy (LB) platform employing low-pass whole genome sequencing (LP-WGS) and targeted sequencing of cell-free (cf) DNA from plasma to detect genome-wide copy number alterations (CNAs) and gene fusions in pediatric solid tumors. A total of 143 plasma samples were analyzed from 19 controls and 73 patients, including 44 bone or soft-tissue sarcomas and 12 renal, 10 germ cell, five hepatic, and two thyroid tumors. cfDNA was isolated from plasma collected at diagnosis, during and after therapy, and/or at relapse. Twenty-six of 37 (70%) patients enrolled at diagnosis without prior therapy (radiation, surgery, or chemotherapy) had circulating tumor DNA (ctDNA), based on the detection of CNAs from LP-WGS, including 18 of 27 (67%) patients with localized disease and eight of 10 (80%) patients with metastatic disease. None of the controls had detectable somatic CNAs. There was a high concordance of CNAs identified by LP-WGS to CNAs detected by chromosomal microarray analysis in the matching tumors. Mutations identified in tumor samples with our next-generation sequencing (NGS) panel, OncoKids®, were also detected by LP-WGS of ctDNA in 14 of 26 plasma samples. Finally, we developed a hybridization-based capture panel to target EWSR1 and FOXO1 fusions from patients with Ewing sarcoma or alveolar rhabdomyosarcoma (ARMS), respectively. Fusions were detected in the plasma from 10 of 12 patients with Ewing sarcoma and in two of two patients with ARMS. Combined, these data demonstrate the clinical applicability of our LB platform to evaluate pediatric patients with a variety of solid tumors.

Potential of Aqueous Humor as a Liquid Biopsy for Uveal Melanoma.

Tumor biopsy can identify prognostic biomarkers for metastatic uveal melanoma (UM), however aqueous humor (AH) liquid biopsy may serve as an adjunct. This study investigated whether the AH of UM eyes has sufficient circulating tumor DNA (ctDNA) to perform genetic analysis. This is a case series of 37 AH samples, taken before or after radiation, and one tumor wash sample, from 12 choroidal and 8 ciliary body (CB) melanoma eyes. AH was analyzed for nucleic acid concentrations. AH DNA and one tumor wash sample underwent shallow whole-genome sequencing followed by Illumina sequencing to detect somatic copy number alterations (SCNAs). Four post-radiation AH underwent targeted sequencing of BAP1 and GNAQ genes. Post-radiation AH had significantly higher DNA and miRNA concentrations than paired pre-radiation samples. Highly recurrent UM SCNAs were identified in 0/11 post-radiation choroidal and 6/8 post-radiation CB AH. SCNAs were highly concordant in a CB post-radiation AH with its matched tumor (r = 0.978). BAP1 or GNAQ variants were detected in 3/4 post-radiation AH samples. AH is a source of ctDNA in UM eyes, particularly in post-radiation CB eyes. For the first time, UM SCNAs and mutations were identified in AH-derived ctDNA. Suggesting that AH can serve as a liquid biopsy for UM.

The spectrum of mitochondrial DNA (mtDNA) mutations in pediatric CNS tumors.

BACKGROUND: We previously established the landscape of mitochondrial DNA (mtDNA) mutations in 23 subtypes of pediatric malignancies, characterized mtDNA mutation profiles among these subtypes, and provided statistically significant evidence for a contributory role of mtDNA mutations to pediatric malignancies. METHODS: To further delineate the spectrum of mtDNA mutations in pediatric central nervous system (CNS) tumors, we analyzed 545 tumor-normal paired whole-genome sequencing datasets from the Children's Brain Tumor Tissue Consortium. RESULTS: Germline mtDNA variants were used to determine the haplogroup, and maternal ancestry, which was not significantly different among tumor types. Among 166 (30.5%) tumors we detected 220 somatic mtDNA mutations, primarily missense mutations (36.8%), as well as 22 loss-of-function mutations. Different pediatric CNS tumor subtypes had distinct mtDNA mutation profiles. The number of mtDNA mutations per tumor ranged from 0.20 (dysembryoplastic neuroepithelial tumor [DNET]) to 0.75 (meningiomas). The average heteroplasmy was 10.7%, ranging from 4.6% in atypical teratoid/rhabdoid tumor (AT/RT) to 26% in diffuse intrinsic pontine glioma. High-grade gliomas had a significant higher number of mtDNA mutations per sample than low-grade gliomas (0.6 vs 0.27) (P = .004), with almost twice as many missense mtDNA mutations per sample (0.24 vs 0.11), and higher average heteroplasmy levels (16% vs 10%). Recurrent mtDNA mutations may represent hotspots which may serve as biologic markers of disease. CONCLUSIONS: Our findings demonstrate varying contributions of mtDNA mutations in different subtypes of CNS tumors. Sequencing the mtDNA genome may ultimately be used to characterize CNS tumors at diagnosis and monitor disease progression.

High Prevalence of SARS-CoV-2 Genetic Variation and D614G Mutation in Pediatric Patients With COVID-19.

BACKGROUND: The full spectrum of the disease phenotype and viral genotype of coronavirus disease 2019 (COVID-19) have yet to be thoroughly explored in children. Here, we analyze the relationships between viral genetic variants and clinical characteristics in children. METHODS: Whole-genome sequencing was performed on respiratory specimens collected for all SARS-CoV-2-positive children (n = 141) between March 13 and June 16, 2020. Viral genetic variations across the SARS-CoV-2 genome were identified and investigated to evaluate genomic correlates of disease severity. RESULTS: Higher viral load was detected in symptomatic patients (P = .0007) and in children <5 years old (P = .0004). Genomic analysis revealed a mean pairwise difference of 10.8 single nucleotide variants (SNVs), and the majority (55.4%) of SNVs led to an amino acid change in the viral proteins. The D614G mutation in the spike protein was present in 99.3% of the isolates. The calculated viral mutational rate of 22.2 substitutions/year contrasts the 13.5 substitutions/year observed in California isolates without the D614G mutation. Phylogenetic clade 20C was associated with severe cases of COVID-19 (odds ratio, 6.95; P = .0467). Epidemiological investigation revealed major representation of 3 of 5 major Nextstrain clades (20A, 20B, and 20C) consistent with multiple introductions of SARS-CoV-2 in Southern California. CONCLUSIONS: Genomic evaluation demonstrated greater than expected genetic diversity, presence of the D614G mutation, increased mutation rate, and evidence of multiple introductions of SARS-CoV-2 into Southern California. Our findings suggest a possible association of phylogenetic clade 20C with severe disease, but small sample size precludes a definitive conclusion. Our study warrants larger and multi-institutional genomic evaluation and has implications for infection control practices.

Increased viral variants in children and young adults with impaired humoral immunity and persistent SARS-CoV-2 infection: A consecutive case series.

BACKGROUND: There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses. METHODS: We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape. FINDINGS: We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape. INTERPRETATION: Our results highlight the potential need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered.

Persistent SARS-CoV-2 infection and increasing viral variants in children and young adults with impaired humoral immunity.

Background: There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses. Methods: We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape. Findings: We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape. Interpretation: Our results highlight the need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered. Funding: The work was partially funded by The Saban Research Institute at Children's Hospital Los Angeles intramural support for COVID-19 Directed Research (X.G. and J.D.B.), the Johns Hopkins Center of Excellence in Influenza Research and Surveillance HHSN272201400007C (A.P.), NIH/NIAID R01AI127877 (S.D.B.), NIH/NIAID R01AI130398 (S.D.B.), NIH 1U54CA260517 (S.D.B.), an endowment to S.D.B. from the Crown Family Foundation, an Early Postdoc.Mobility Fellowship Stipend to O.F.W. from the Swiss National Science Foundation (SNSF), and a Coulter COVID-19 Rapid Response Award to S.D.B. L.G. is a SHARE Research Fellow in Pediatric Hematology-Oncology.

Early pandemic molecular diversity of SARS-CoV-2 in children.

Background: In the US, community circulation of the SARS-CoV-2 virus likely began in February 2020 after mostly travel-related cases. Children's Hospital of Philadelphia began testing on 3/9/2020 for pediatric and adult patients, and for all admitted patients on 4/1/2020, allowing an early glimpse into the local molecular epidemiology of the virus. Methods: We obtained 169 SARS-CoV-2 samples (83 from patients <21 years old) from March through May and produced whole genome sequences. We used genotyping tools to track variants over time and to test for possible genotype associated clinical presentations and outcomes in children. Results: Our analysis uncovered 13 major lineages that changed in relative abundance as cases peaked in mid-April in Philadelphia. We detected at least 6 introductions of distinct viral variants into the population. As a group, children had more diverse virus genotypes than the adults tested. No strong differences in clinical variables were associated with genotypes. Conclusions: Whole genome analysis revealed unexpected diversity, and distinct circulating viral variants within the initial peak of cases in Philadelphia. Most introductions appeared to be local from nearby states. Although limited by sample size, we found no evidence that different genotypes had different clinical impacts in children in this study.

Resistance to autosomal dominant Alzheimer's disease in an APOE3 Christchurch homozygote: a case report.

We identified a PSEN1 (presenilin 1) mutation carrier from the world's largest autosomal dominant Alzheimer's disease kindred, who did not develop mild cognitive impairment until her seventies, three decades after the expected age of clinical onset. The individual had two copies of the APOE3 Christchurch (R136S) mutation, unusually high brain amyloid levels and limited tau and neurodegenerative measurements. Our findings have implications for the role of APOE in the pathogenesis, treatment and prevention of Alzheimer's disease.

A semiautomated whole-exome sequencing workflow leads to increased diagnostic yield and identification of novel candidate variants.

Advancing the clinical utility of whole-exome sequencing (WES) for patients with suspected genetic disorders is largely driven by bioinformatics approaches that streamline data processing and analysis. Herein, we describe our experience with implementing a semiautomated and phenotype-driven WES diagnostic workflow, incorporating both the DRAGEN pipeline and the Exomiser variant prioritization tool, at an academic children's hospital with an ethnically diverse pediatric patient population. We achieved a 41% molecular diagnostic rate for 66 duo-, quad-, or trio-WES cases, and 28% for 40 singleton-WES cases. Preliminary results were returned to ordering physicians within 1 wk for 12 of 38 (32%) probands with positive findings, which were instrumental in guiding the appropriate clinical management for a variety of patients, especially in critical care settings. The semiautomated and streamlined WES workflow also enabled us to identify novel variants in candidate disease genes in patients with developmental delay and autism and immune disorders and cancer, including ANK2, BPTF, BCL11A, FOXN1, PLAA, ATRX, DNAJC21, and RAD50 Together, we demonstrated the implementation of a streamlined WES workflow that was successfully applied for both clinical and research purposes.

OncoKids: A Comprehensive Next-Generation Sequencing Panel for Pediatric Malignancies.

The OncoKids panel is an amplification-based next-generation sequencing assay designed to detect diagnostic, prognostic, and therapeutic markers across the spectrum of pediatric malignancies, including leukemias, sarcomas, brain tumors, and embryonal tumors. This panel uses low input amounts of DNA (20 ng) and RNA (20 ng) and is compatible with formalin-fixed, paraffin-embedded and frozen tissue, bone marrow, and peripheral blood. The DNA content of this panel covers the full coding regions of 44 cancer predisposition loci, tumor suppressor genes, and oncogenes; hotspots for mutations in 82 genes; and amplification events in 24 genes. The RNA content includes 1421 targeted gene fusions. We describe the validation of this panel by using a large cohort of 192 unique clinical samples that included a wide range of tumor types and alterations. Robust performance was observed for analytical sensitivity, reproducibility, and limit of detection studies. The results from this study support the use of OncoKids for routine clinical testing of a wide variety of pediatric malignancies.

Inherited germline ATRX mutation in two brothers with ATR-X syndrome and osteosarcoma.

We report a family in which two brothers had an undiagnosed genetic disorder comprised of dysmorphic features, microcephaly, severe intellectual disability (non-verbal), mild anemia, and cryptorchidism. Both developed osteosarcoma. Trio exome sequencing (using blood samples from the younger brother and both parents) was performed and a nonsense NM_000489.4:c.7156C>T (p.Arg2386*) mutation in the ATRX gene was identified in the proband (hemizygous) and in the mother's peripheral blood DNA (heterozygous). The mother is healthy, does not exhibit any clinical manifestations of ATR-X syndrome and there was no family history of cancer. The same hemizygous pathogenic variant was confirmed in the affected older brother's skin tissue by subsequent Sanger sequencing. Chromosomal microarray studies of both brothers' osteosarcomas revealed complex copy number alterations consistent with the clinical diagnosis of osteosarcoma. Recently, somatic mutations in the ATRX gene have been observed as recurrent alterations in both osteosarcoma and brain tumors. However, it is unclear if there is any association between osteosarcoma and germline ATRX mutations, specifically in patients with constitutional ATR-X syndrome. This is the first report of osteosarcoma diagnosed in two males with ATR-X syndrome, suggesting a potential increased risk for cancer in patients with this disorder.

Generation and characterization of tamoxifen-inducible Pax9-CreER knock-in mice using CrispR/Cas9.

Pax9 encodes a paired-box homeodomain (Pax) transcription factor and is critical for the development of multiple organs. Using CrispR/Cas9-mediated homologous directed repair (HDR), we generated a new Pax9-CreER knock-in mouse line in which the CreER(T2) fusion protein is produced after synthesis of endogenous Pax9 protein. We found that tdTomato reporter expression in Pax9-CreER;tdTomato reporter mice is detectable in a similar pattern to the endogenous Pax9 expression, faithfully recapitulating the Pax9 expression domains throughout the embryo and in the adult mouse. At early embryonic stages, the tdTomato reporter is expressed first in the pharyngeal pouch region and later in the craniofacial mesenchyme, somites, limbs, and lingual papillae in the adult tongue. These results demonstrate that this new Pax9-CreER knock-in mouse line can be used for lineage tracing and genetic targeting of Pax9-expressing cells and their progeny in a temporally and spatially controlled manner during development and organogenesis.

Comprehensive Molecular Characterization of Papillary Renal-Cell Carcinoma.

BACKGROUND: Papillary renal-cell carcinoma, which accounts for 15 to 20% of renal-cell carcinomas, is a heterogeneous disease that consists of various types of renal cancer, including tumors with indolent, multifocal presentation and solitary tumors with an aggressive, highly lethal phenotype. Little is known about the genetic basis of sporadic papillary renal-cell carcinoma, and no effective forms of therapy for advanced disease exist. METHODS: We performed comprehensive molecular characterization of 161 primary papillary renal-cell carcinomas, using whole-exome sequencing, copy-number analysis, messenger RNA and microRNA sequencing, DNA-methylation analysis, and proteomic analysis. RESULTS: Type 1 and type 2 papillary renal-cell carcinomas were shown to be different types of renal cancer characterized by specific genetic alterations, with type 2 further classified into three individual subgroups on the basis of molecular differences associated with patient survival. Type 1 tumors were associated with MET alterations, whereas type 2 tumors were characterized by CDKN2A silencing, SETD2 mutations, TFE3 fusions, and increased expression of the NRF2-antioxidant response element (ARE) pathway. A CpG island methylator phenotype (CIMP) was observed in a distinct subgroup of type 2 papillary renal-cell carcinomas that was characterized by poor survival and mutation of the gene encoding fumarate hydratase (FH). CONCLUSIONS: Type 1 and type 2 papillary renal-cell carcinomas were shown to be clinically and biologically distinct. Alterations in the MET pathway were associated with type 1, and activation of the NRF2-ARE pathway was associated with type 2; CDKN2A loss and CIMP in type 2 conveyed a poor prognosis. Furthermore, type 2 papillary renal-cell carcinoma consisted of at least three subtypes based on molecular and phenotypic features. (Funded by the National Institutes of Health.).